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The subject discussed in this text is the oscillatory behavior of solutions of holomorphic linear 
time-dependent differential systems of equations. In [1] a bound is provided for the number of 
times a trajectory of a polynomial vector field crosses a given algebraic hypersurface. This implies 
a bound also for linear time dependent vector fields with polynomial dependence on time. A dif- 
ferent approach is taken here, giving better bound for some special cases. 

This text is a preliminary report. 

Besides of having an intrinsic interest, a related problem arises in estimating the number of zeros 
which the so called complete Abelian integrals may have, the latter problem being itself connected 
to the problem of estimating the number of limit cycles which can be born from polynomial per- 
turbations of Hamiltonian (polynomial) vector fields on the plane (see, for example, [4]). 

Remark. The bounds are given here for a real segment in the domain where the coefficients of 
the equations are holomorphic and real on the real axis, but they can be given in principle for any 
compact domain contained in the domain where the coefficients are holomorphic ([3]). We also 
remark that the results on which the bounds appearing here are based, are not the best available 
(see discussion in [3]). 

Let us be given a linear differential system of equations 

x = A(t)x 

with coefficients holomorphic in some simply connected complex domain. Then the solution will 
be a vector valued holomorphic function in this domain. The problem in the large is to estimate 
the number of zeros the first component of the solution may have in that domain. 

It is known that if the coefficients of the system belong to the class of polynomials of a given 
degree, a constructive bound on the number of zeros on some compact subset of the domain may 
be given which involves just the degree and the maximum modulus of the coefficients [1]. 



In the following, another way of obtaining bounds on the number of zeros of the solutions of 
such systems is discussed, which is the consequence of a simple but somewhat surprising fact. 

In [1] , one of the aims is to bound the complex oscillation of solutions of polynomial vector fields; the 
result for linear differential systems is obtained as a consequence. One writes the system as a poly- 
nomial differential system of equations with the coefficients and time treated as extra variables, 
the polynomial coefficients being integral. Then, by subsequent differentiating, one obtains the 
derivatives of the first variable written as polynomials with integer coefficients in all the variables 
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of the system. It turns out that one can effectively bound the moment when the ideal generated 
by the latter polynomials stabilizes. This last derivative is then written as a linear combination of 
the former derivatives over the ring of polynomials in the variables of the system, which is a linear 
scalar differential equation for the first coordinate. Since the coefficients of the polynomials are 
integers, one uses a form of effective Nullstellensatz to bound the coefficients of this equation (on 
compact regions in the variable space) . Since for linear scalar differential equations a bound on the 
number of zeros exists in terms of the maximum of the coefficients in the region ( [2] ) , a bound on 
the number of zeros is obtained. 

One notes that the obtained scalar differential equation has a space of solutions of dimension 
much larger than the dimension of the original system, and thus has many solutions not corre- 
sponding (at least directly) to the solutions of the system. 

In what we do here, we keep the integrality of the polynomial coefficients to be able to use the 
effective Nullstellensatz. As in [1], we can parameterize in this way all linear differential systems 
with coefficients being polynomials of a fixed degree. In the beginning, however, we assume only 
holomorphic dependence on the time and parameters so that not to obscure the main point. 

So let there be given a linear differential system of order n x = A(t,p)x = (aij)(t,p)(x), where 
p stands for parameters. We assume the dependence of the system matrix on t,p to be holomorphic 
in some region U. By W we denote the projection of the domain on the parameter space. We 
describe a procedure to obtain a scalar differential equation for the first system variable. Differen- 
tiating and substituting the relation we have for the first variable, we get: 

n n n n n 

x" = ^2 a u xi + aux'i = ^2 a 'u x i + a u ^2 ai i x i = X^ ' 1 * + a ij a ji) x i- 

i=l i=l i=l i=l j=l 

(i) 

Similarly we get the expressions for x\ for any i. We write 

x\ = y cij Xj = a y 'x. 

There is a minimal k, 1 < k < n, such that a^ \a^\ ..,a^ are linearly dependent. We then have 
a unique decomposition 

a^=A k . ia ^ + ... + AW\ 

A,i £ M.{U), implying that the following equation must be satisfied by any solution for the x\ 
variable of the initial system 

xP -A fc _ixf _1) -...-A QXl = 0. 

In the sequel we call this differential equation simply 'the derived equation'. Notice that the 
determinant of a suitable fe-minor of the matrix 

(aWaW.-.a**- 1 )) 

is not identically zero, and that multiplying by it the equation, we get an equation with holomor- 
phic coefficients. We thus call it the leading coefficient of the derived equation (corresponding to 
the minor chosen; of course in the case k=n there is only one such choice). 
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Claim 1. There exists an analytic set G such that for any p G W — G any solution of the above 
equation is a solution for the x\ variable of the initial system. 

Proof. The leading coefficient of the derived equation is a holomorphic function in p, t, not 
identically zero. Thus the values of p for which the leading coefficient vanishes identically in t, 
form an analytic subset of W, G. For p E W — G for all but a discrete set of values of t we have 
rank (a^\t, e) = k. 

Without limiting generality, we assume U to be a polydisc. For any fixed po € W — G, there 
is a natural mapping from the space of solutions of the system to the space of solutions of the 
derived equation, which we can identify with the linear mapping 

xi(t ,po),...,x n (t ,po) i * x 1 (t ,p ),...,x[ k) (t ,p ) 
where to is such that the leading coefficient does not vanish on to,po, given by 

Xi\to,Po) = (4 J) (io,Po))(^j(io,Po))- 

Since the leading coefficient does not vanish at to,po, rank (a^\t,e)) = k, and the mapping is 
surjective. Hence Vpo E W — G every solution of the derived equation comes from a solution of the 
system. 

Remark: if the dependence is algebraic rather than holomorphic, the set G is algebraic. 
Let us express the determinant of the i-th. A;-minor, the order of the derived equation being k, as 
ho{p) + bn{p)t + ... + bik t t k \ We call the ideal generated by all bij(p) the degeneracy ideal, and the 
corresponding variety the degeneracy variety. 

We now make a small digression, and explain a basic difficulty which might occur in estimat- 
ing the number of zeros of the first component of a solution of the system through the derived 
above equation for this component. In [2] one finds the following theorem, derived from growth 
estimates on the solutions: 

Theorem. Let 

a(t)y (n) + an-i^y*"- 1 ) + ... + a (t)y = 

have coefficients holomorphic in a region U containing a real segment K . Suppose that the modulus 
of all ai(t) is bounded in U by A, and that 

maxt^K\a(t)\ > a. 

Then the number of zeros on K of any holomorphic in U solution of the equation, real on the real 
axis, is not greater than 

{A/a + nY, 

where fi depends only on the geometry of U and K . 

The main obstacle in applying this theorem to our case lies in the fact that, as parameters vary, 
maxt<=K\a(t,p)\ (which depends on the parameters) a priori may become arbitrarily small - the 
situation of singular perturbations - and the estimates explode. 
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Our next goal is to show that the scalar differential equation we derived for x\ cannot be in 
fact singularly perturbed, at least when one parameter is involved. This is somewhat unexpected, 
and we explain why. Consider the following very simple e-dependent system: 

x = x + ey y = x + y 

0, the derived equation, which is of second order, will simply 
x, implying that the derived equation is singularly perturbed. 
(1 + e)x + 2ey, and the derived equation is 

x - 2x + (1 -e)x = 0. 

This equation is not singularly perturbed. It has the property that for the exceptional value of the 
parameter, an additional solution appears which does not come from the original system, namely 
te t . To explain why the singular perturbations cannot occur also in the general case, we proceed as 
follows. Until further notice the dimension of the parameter space is one, and the parameter will 
be denoted by e. 

Definition. We call a meromorphic function f(t,e) singularly perturbed at eo, if when expressed 
as a quotient of two holomorphic functions f(t,e) = u(t,e)/v(t,e), u(t,e) is divisible by a higher 
or equal power of (e — £o) than the one that divides v (t, e). 

Definition. A (scalar) differential equation with meromorphic coefficients depending on one pa- 
rameter, normalized so that the leading coefficient is equal to 1, will be called singularly perturbed 
if at least one of its coefficients is singularly perturbed at some point in the parameter space. 

We make the following obvious claim. 

Lemma 1. Let us be given a differential equation with t, e-dependent coefficients meromorphic 
in a domain U , singularly perturbed at Eq. Then there is a point to, such that in the neighborhood 
of(to,£o) £ U, the equation takes the form 

(e - £ )V n) + a n _!(t, e)y^ + ... + a (t, e)y = 0, 

s > 0, with coefficients holomorphic in the neighborhood of(to,Eo), and at least one of them is not 
identically zero when e = Eq. 

Theorem 1. Let a differential equation of order n with coefficients which depend meromorphi- 
cally on t,e in some polydisc U x V, be singularly perturbed for at least one parameter value Sq, 
and let it have I solutions yi(t, e), yi{t, e) holomorphic in t and e in that polydisc. Then for each 
fixed parameter value i from V , 

dime span c (yi(t,i), ...,yi(t,i)) < n-1. 

Proof. Using Lemma 1, we write the equation in the form given there in some smaller polydisc, 
and assume to = 0, Sq = 0: 

e s y {n) + on-ift e)y^ l) + ... + a (t, e)y = 0, 

with afc(0, 0) 7^ 0. We claim that if two holomorphic in t, e solutions are such that they have equal 
y( n ~ 1 )(0,e), ...,y( fc+1 )(0, e),y( A: ~ 1 )(0,£), ...,y(0,e), they are in fact identical. Indeed, otherwise we 
get that there exists a nonzero holomorphic in t,e solution y(t,e) such that 

y^^e) = ... = y( k+1 H0,e) = y^^e) = ... = y(0,e) = 0. 



One might expect that when e = 
degenerate into the first order x = 
However, x = x + ey + e(x + y) = 
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Since y(t,e) is not identically zero, y' ^(0, e) cannot be identically zero. Thus y (fc) (0,e) = e h g(e), 
h > 0, g(0) / 0. Assume h > 0. y(t,0) is the solution of the differential equation with e set to 
zero. Since a*; (0,0) / and the initial conditions at e = are zero (h > 0), we have y(t, 0) = 
identically. This can only happen if y(t,e) = ey±(t, e), yi(t, e) being also a holomorphic solution of 
the system with the same zero initial conditions except for y^(0,e) = e h ~ 1 g(e). We continue in 
the same way, until we get a holomorphic solution yh(t, e) = y(t,e)/e h with the initial conditions 

y^X^e) = ... = y (k+1 H0,e) = y^^e) = ... = y(0,e) = 0, 

and y( k \o,e) = g(e), g(0) / 0. However substituting yh(t,e) into the original equation and setting 
t = 0, we get: 

e s y^(0,e)+a k (0,e)y^(0,e)=0, 
which is a contradiction (ctfc(0, 0), y^ (0, 0) / 0). 

Thus in the smaller polydisc the holomorphic in t, e solutions can be identified with sets of holomor- 
phically dependent on e sets of initial conditions y( n_1 ^(0, e), ?/ fc+1 ) (0, e), y^ -1 ) (0, e), ...,y(0, e), 
omitting y( k \o, e). Now we prove the theorem. Suppose we have holomorphic in U x V solutions 
of our singularly perturbed equation yi(t,e), ..,yi(t,e). Over the field meromorphic in V functions 
there can be at most n — 1 linearly independent such solutions since the latter are identified with 
0{V) n ~ l . So we choose a maximal set S of linearly independent over Ai(V) solutions and express 
all the other as their linear combination over A4(V). It follows that all the other solutions lie in the 
space spanned by the solutions from S for all but a discrete set of e-s. Thus for almost all values of 
e, the space of holomorphic in t, e solutions is n — 1 dimensional, which implies that it is at most 
n — 1 dimensional for all values of e by the semicontinuity of the dimension. 

Theorem 2. Suppose that we are given a linear differential system with coefficients depending 
on one parameter, holomorphic both in t and e in a polydisc U x V. Then the derived scalar dif- 
ferential equation for the variable x\ cannot be singularly perturbed. 

Proof. Returning to notations of proof of Claim 1, we take a point (io> £ o) £U — Z, m which the 
leading coefficient of the derived equation is not zero, and rank {a^\t,e) = k. 

For each m, there is a holomorphic in U solution of the system defined by the initial conditions 
(constant in e): 

xi(t ,e) = ... x m (t ,e) = 1 ... x n (t ,e) = 0. 

We denote the 1-st component of this solution by x™(t,e). By discussion from the proof of Claim 
1, since x\(t, eq), x"(t, £q) span the space of solutions of the derived equation at eo, and since at 
this (generic) point the order of the derived equation is k, we have 

dime spanc x\(t,£ ), ...,Xi(t,£ ) = k. 

However, if the derived equation is singularly perturbed even at only one point E\ G V, from 
Theorem 1 it would follow that for every eq £ V 

dime spanc x\(t,£ ), ...,x i (t,e ) < k — 1, 
giving a contradiction. Thus the derived equation cannot be singularly perturbed. 
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We explain now our goal in what follows. As it was already mentioned, we can parameterize 
the problem of general linear system of order n with polynomial coefficients of degree d as: 

n d 

Xi = y ^ ^ o,ijkt Xj i = 1, .., n 

j=l fc=0 

the parameters being a^, thus making the system matrix to be polynomial in the time and pa- 
rameters, with integer coefficients. Ideally we would like to bound the number of zeros which the 
first variable of this system might have, but the results available here are still not permitting such 
general bounds. Instead, we restrict ourselves to the case when the system matrix is polynomial in 
time and one additional parameter, then making several remarks about the general linear differen- 
tial system with polynomial coefficients as above. 

So let us have a linear differential system of order n, with coefficients being polynomials in t, e of de- 
gree d, with integer coefficients of maximum modulus M. Let the order of the derived equation be k. 

Lemma 2. Let a nonsingularly perturbed rational function be written as 

a (e) + ai (e)t + ... + a k {e)t k 
b {e)+b 1 {e)t + ... + b l (e)t l ' 

Then each aj belongs to the ideal generated by bo , b±, .., b\. 

Proof, aj must vanish at points where all b±,..,bi vanish (otherwise the fraction would be sin- 
gularly perturbed). Write bi(e) = Ci{e)b(e), where b{e) is the common factor. Then Cj(e) have no 
common zero and therefore, since we are working over C, 1 = hi(e)c\{e) + .. + hi(e)ci(e), for some 
polynomials hi,.., hi. But aj must be divisible by b(e): otherwise we get singular perturbations. 
Thus: 

oj(e) = d{e)b{e) = d{e)b{e){h 1 {e)c 1 {e) + .. + hi(e)ci{e)) = d(e)hi(e)bi(e) + .. + d(e)hi(e)bi(e) . 
We derive now some computational estimates we will need. 

Lemma 3. Suppose we are given a non- singularly perturbed rational function oft,e, 

q (g) + ai{e)t + ... + a k {e)t k 
b (e) + b 1 (e)t + ... + b l (e)t l ' 

with degrees of the numerator and the denominator being at most d, and the coefficients being 
integers of maximum modulus M (we mean by this possibly complex integers). Then one may write: 

ao(£)+a 1 (e)t + ... + a k (e)t k _ giU j Ej=i Mf)Mf) 
6o(e) + h(e)t + ... + h{e)t l b (e) + h(e)t + ... + h{e)t l ' 

with hij(e) being polynomials of degree at most 2d — 1 with integer coefficients of maximum modulus 
{eM) cS for some constant C . 

Proof. From Lemma 2, a« = ^2hijbj. By [5], since Vi deg(ai) < d, Vj deg(bj) < d, we have 
that aj = h'ijbj-, degQiij) <2d—l. Observe that this gives us an integer coefficient linear system 
of dimension at most 2d{d+ 1), maximum modulus of the coefficients not greater than M. By the 
integrality trick this gives us the estimate (2d(d + \))\M 2d ^ d+l ^ < (eM) Cd3 on the coefficients of 
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the polynomials hij. 

Lemma 4. In the setting of Lemma 3, suppose that \e\ is bounded by E. Then the above fraction 
can always be written as 

Eti *«* 

with a < (eM) CdS (E^- 1 + 1) for some constant C , with some di equal to 1 and all \d{\ being not 
greater than 1. 

Proof. For each e in the disc, take the bi having the maximum modulus, and divide by it both 
the numerator and the denominator: use Lemma 3 (the constant here is different from the one in 
Lemma 3). 

Lemma 5. Suppose we are given a non- singularly perturbed differential equation of order k 

(k) , Qfc-i(*,e) (fc-i) , , ^M,. _ n 
V + b(t,e) V + - + b(t,e) V ' 

all cii(t, e) and b(t, e) being polynomials of degree at most d, having integer coefficients of maximum 
modulus M . 

Let [-R/2, R/2] be a segment on the real axis containing [—1,1]. Fix any real e of modulus less 
than E. Let f be a holomorphic solution of the equation for this parameter value, real on the real 
axis. Then the number of zeros f can have on {—R/2, R/2], is bounded by 

({Me) Cd \E 2d + l)R d+1 + ky , 

with a and C being some constants. 

Proof. By Lemma 4, we can write the equation in the form 

(t s + J2 diWv™ + (cfc-i o + •• + c fc -i d t d )y^ + ... + (coo + - + c 0d t d )y = 0. 

with |cjj| being bounded by the bound from Lemma 2, and all < 1. If we could find now some 
point in [— R/2, R/2] where t s + J2i^ s ^ 1 nas some value bounded away from zero, we are done. 
This because then we may just apply the Theorem cited from [2]. 
So we write 

+ > \t s + d s _ 1 t s - 1 + ... + d \ - \t s+1 (d s+1 + d s+2 t + ..)\ 

and notice that by Cartan lemma ([2]), one may delete disks of total diameter h from C, such that 
on the rest of C the modulus of the monic polynomial t s + dg^it 3 ^ 1 + ... + do is > (/i/4e) s . Thus 
for each h > there exists a point th of modulus at most h/2, at which the value of the polynomial 
is > (/i/4e) s . But then: 

\t s h + J2 d ^h\ Z \t s h + d s - 1 t s - 1 + ... + d \-\t s h +1 (d s+1 + d s+2 t h + ..)\ > (h/4e) s -h s+1 (d-s). 
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Taking for example h = l/(2(4e) s (d — s)) we get that there exists a point on the unit segment at 
which \t s + Yli^s dit l \ is bounded from below by 



> 



(4e) s + s 2 s + 1 (d- s) s ~ (4e) s 2 + s 2 s + 1 d s 

As s can be from to d — 1 (it can be shown that the case s = d is also covered by the obtained 
estimate), we are led to the conclusion that there will be a point in [-R/2, R/2] where 



\t s +j2 d i f \ > u e) <p+d 2 d+i d d- 



In turn, | ^j =0 CijP\ cannot, when e varies over disk of radius E and t over Dr, be larger than 

d 

Y;R j ■ {eM) Cd \E 2d - 1 + 1) < {eM) Cd \E 2d - 1 + l)R d+1 . 

j=0 

By Theorem cited from [2], one deduces that for any given real |e| < E, any holomorphic in t 
solution of the given differential equation, real on the real axis, may have not more than 

/ {eM) Cd3 (E 2d ~ l + l)R d+1 \ a 
\ \/{Ae) d2 + d 2 d +H d + ) 

zeros on [— R/2, R/2], a not depending on R because of properties of a{K,U) from the Theorem. 
Simplifying, one gets the bound we stated. 

To estimate the number of zeros a solution (more precisely, its x\ component) of the initial linear 
system 

x = A(t, e)x 

might have, we now have only to bound the maximum modulus of the (integer) coefficients of the 
derived equation. We will do the computation for an arbitrary number of parameters; to emphasize 
this we shall write p instead of e. To remind the notations, we denote the dimension of the system by 
n, the degree of the polynomial coefficients by d, the maximum modulus of the coefficients of those 
polynomials by M, the number of parameters by q, and the order of the derived equation for x\ by k. 

Lemma 6. The following recurrent formulas hold for oS- l \ defined in the beginning of this text: 

a W =a^ +A(t,p) T a^ 1 \ 

Proof. Immediate. 



The following is a computation: 



Lemma 7. Let Pi, Pi be two polynomials of degrees d\,d2 in s variables, with the modulus of 
their coefficients being not larger than M\,M2, respectively. Then the degree of P1P2 is d\ +di and 
the coefficients are not larger in absolute value than (1 + min(d\, d2)) s M1M2. 

Lemma 8. The degrees and the maximum modulus of coefficients of the entries of are di 
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and n l (d + (d + l) q+l M) 1 , respectively. 



Proof. By writing recursive relations. 



We now compute the maximum modulus and the degrees for the coefficients of the derived equation. 
Lemma 9. Let 

{k) a k -i(t,p) (k -D (*o(t,p) = 

U P(t,p) U P(t,p) U 

be the equation derived for x\ from the original system. Then the degrees of both a, and (3 are 
bounded by k(k + l)d/2, and the maximum modulus for the (integer) coefficients of those polyno- 
mials is not larger than (eM(d + l)) c ^ +1 ^ n for some constant C. 



Proof. Again by recursive relations. 



Theorem 2. Let the initial linear differential system have as coefficients polynomials in t,e (e one 
dimensional) with integer coefficients of maximum modulus M, and let [-R/2, R/2] be a segment 
on the real axis containing [—1,1]. Fix any real e of modulus less than E. Let f be a holomorphic 
solution of the equation for this parameter value, real on the real axis. Then the number of zeros f 
can have on [— R/2, R/2], is bounded by 

((Me) Cn9rf4 (£ n(n+1)d + i)R^(n+i)d/2+i + n y ^ 
for some constants C and a. 

Proof. Just substitute the bounds stated in Lemma 9 with q = 1 into Lemma 5 and simplify 
the expression. 

Of course, the concrete form of the bound given above is irrelevant; what is important is the 
rate of its growth, which is simple exponential in n, d compared with the tower of exponents which 
can be derived from [1] in this situation. 

We now make some remarks on the general case when many parameters are allowed. First, it is 
clear that the absence of singular perturbations in the case of dependence on one parameter implies 
absence of singular perturbations along any given holomorphic curve in the parameter space, and 
the following is true. 

Lemma 2a. Let a rational function dependent on several parameters be written as 

ao{p) + ai(p)t + ••■ + (ik(p)t k 
b {p) + b 1 {p)t + ... + b l (p)t l ' 

and let it be such that it is not singularly perturbed along any line in the parameter space. Then 
each a,j vanishes on the variety where all bo,b±,..,bi vanish. 

We cannot in general have the conclusion of Lemma 2 as the example of the nonsingularly perturbed 
fraction 

ab 
a 2 + b 2 t 
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shows. That is, in the derived equation for system with polynomial coefficients depending on many 
parameters we can only guarantee that the coefficients in the numerator vanish on the degeneracy 
variety, but not that they in fact belong to the degeneracy ideal. However, if by some additional 
reasoning for, say, the general linear system with polynomial coefficients 

n d 

%i = ^ ^ ^ ^ O'ijkt Xj i = 1, .., 77., 
j=l k=0 

it will turn out to be true, then the bound, though drastically grows because of the double expo- 
nential in n 2 (d + 1) bounds in the effective Nullstellensatz (see [5]), has the asymptotic form (in 
terms of Theorem 2) 

e e^ d \ ER fn^ 



where P is a polynomial. This is still a considerable improvement compared to [1] (of course if the 
assumption is valid). 



References. 

1. Novikov, D.; Yakovenko, S. Trajectories of polynomial vector fields and ascending chains of 
polynomial ideals. Ann. Inst. Fourier (Grenoble) 49 (1999), no. 2, 563-609. 

2. IPyashenko, Y.; Yakovenko, S. Counting real zeros of analytic functions satisfying linear or- 
dinary differential equations. J. Differential Equations 126 (1996), no. 1, 87-105. 

3. Yakovenko S. On functions and curves defined by ordinary differential equations , to appear 
in Proceedings of the Arnoldfest (Ed. by Bierstone, Khesin, Khovanskii, Marsden), Fields Institute 
Communications, 1998. 

4. Il'yashenko, Y.; Yakovenko, S. Double exponential estimate for the number of zeros of complete 
abelian integrals and rational envelopes of linear ordinary differential equations with an irreducible 
monodromy group. Invent. Math. 121 (1995), no. 3, 613-650. 

5. Shiffman, B. Degree bounds for the division problem in polynomial ideals. Michigan Math. 
J. 36 (1989), no. 2, 163-171. 



10 



